Validity of semiclassical gravity in the stochastic gravity approach 
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In semiclassical gravity the back-reaction of the classical gravitational field interacting with quan- 
tum matter fields is described by the semiclassical Einstein equations. A criterion for the validity of 
semiclassical gravity based on the stability of the solutions of the semiclassical Einstein equations 
with respect to quantum metric perturbations is discussed. The two-point quantum correlation 
functions for the metric perturbations can be described by the Einstein-Langevin equation obtained 
in the framework of stochastic gravity. These correlation functions agree, to leading order in the 
large N limit, with the quantum correlation functions of the theory of gravity interacting with N 
matter fields. The Einstein-Langevin equations exhibit runaway solutions and methods to deal with 
these solutions are discussed. The validity criterion is used to show that flat spacetime as a solution 
of semiclassical gravity is stable and, consequently, a description based on semiclassical gravity is a 
valid approximation in that case. 
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I. INTRODUCTION 

Semiclassical gravity describes the interaction of the 
gravitational field as a classical field with quantum mat- 
ter fields. For a free quantum field this theory is robust 
in the sense that it is consistent and fairly well under- 
stood P, • The gravitational field is described by the 
semiclassical Einstein equation which has as a source the 
expectation value in some quantum state of the matter 
stress tensor operator. The semiclassical theory is in 
some sense unique as a theory where the gravitational 
field is classical. In fact, a classical gravitational field in- 
teracts with other fields through their stress tensors, and 
the only reasonable c-number stress tensor that one may 
construct [3|, |^, |5j with the stress tensor operator of a 
quantum field is its expectation value in some quantum 
state. However, the scope and limits of the theory are not 
so well understood because we still lack a fully well under- 
stood quantum theory of gravity. It is assumed that the 
semiclassical theory should break down at Planck scales, 
which is when simple order of magnitude estimates sug- 
gest that the quantum effects of gravity cannot be ig- 
nored: the gravitational energy of a quantum fluctua- 
tion of energy in a Planck size region, determined by the 
Heisenberg uncertainty principle, is of the same order of 
magnitude as the energy of the fluctuation itself. 

From the semiclassical Einstein equations it seems also 
clear that the semiclassical theory should break down 
when the quantum fluctuations of the stress tensor are 
large. Ford @ was among the first to have emphasized 
the importance of these quantum fluctuations. It is less 
clear, however, how to quantify the size of these fluc- 
tuations. Thus, Kuo and Ford |jj used the variance of 
the fluctuations of the stress tensor operator compared 
to the mean value as a measure of the validity of semi- 
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classical gravity. As pointed out by Hu and Phillips 
0,0 such a criterion should be refined by considering the 
back reaction of those fluctuations on the metric. Ford 
and collaborators also noticed that the metric fluctua- 
tions associated to the matter fluctuations can be mean- 
ingfully classified as "active" [13, HH E3 and "passive" 

SlIilQIIi. 

A different approach to the validity of semiclassical 
gravity was taken by Horowitz 0, 0] who studied the 
stability of a semiclassical solution with respect to lin- 
ear metric perturbations. In the case of a free quan- 
tum matter field in its Minkowski vacuum state, flat 
spacetime is a solution of semiclassical gravity. The 
equations describing those metric perturbations involve 
higher order derivatives, and Horowitz found unstable 
"runaway" solutions that grow exponentially with char- 
acteristic timescales comparable to the Planck time; see 
also the analysis by Jordan [l^. Later, Simon [l9l. |20|. 
argued that those unstable solutions lie beyond the ex- 
pected domain of validity of the theory and emphasized 
that only those solutions which resulted from truncat- 
ing perturbative expansions in terms of the square of the 
Planck length are physically acceptable |19t |2(j ■ Further 
discussion was provided by Flanagan and Wald |2l| , who 
advocated the use of an "order reduction" prescription 
first introduced by Parker and Simon [22| • More recently 
Anderson, Molina-Paris and Mottola have taken up the 
issue of the validity of semiclassical gravity (2^ again. 
Their starting point is the fact that the semiclassical 
Einstein equation will fail to provide a valid description 
of the dynamics of the mean spacetime geometry when- 
ever the higher order radiative corrections to the effective 
action, involving loops of gravitons or internal graviton 
propagators, become important. Next, they argue quali- 
tatively that such higher order radiative corrections can- 
not be neglected if the metric fluctuations grow without 
bound. Finally, they propose a criterion to characterize 
the growth of the metric fluctuations, and hence the va- 
lidity of semiclassical gravity, based on the stability of the 
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solutions of the linearized semiclassical equation. Follow- 
ing these approaches the Minkowski metric is shown to 
be a stable solution of semiclassical gravity with respect 
to small metric perturbations. 

As emphasized in Ref. 23] the above criteria may be 
understood as criteria within semiclassical gravity itself. 
It is certainly true that stability is a necessary condi- 
tion for the validity of a semiclassical solution, but one 
may also look for criteria within extensions of semiclas- 
sical gravity. In the absence of a quantum theory of 
gravity such criteria may be found in more modest ex- 
tensions. Thus, Ford considered graviton production 
in linearized quantum gravity and compared the results 
with the production of gravitational waves in semiclas- 
sical gravity. Ashtekar [24| and Beetle [25| found large 
quantum gravity effects in three-dimensional quantum 
gravity models. In a recent paper [2ij (see also Ref. 
we advocate for a criteria within the stochastic gravity 
approach. Stochastic semiclassical gravity extends semi- 
classical gravity by incorporating the quantum stress ten- 
sor fluctuations of the matter fields; see Refs. (2^, 29] for 
reviews. 

It turns out that this validity criteria is equivalent to 
the validity criteria that one might advocate within the 
large N expansion, that is the theory describing the inter- 
action of the gravitational field with N identical matter 
fields. In the leading order, namely the limit in which 
N goes to infinity and the gravitational constant is ap- 
propriately rescaled, the theory reproduces semiclassical 
gravity. Thus, a natural extension of semiclassical grav- 
ity is provided by the next to leading order. It turns out 
that the symmetrized two-point quantum correlations of 
the metric perturbations in the large N expansion are 
equivalent to the two-point stochastic metric fluctuations 
predicted by stochastic gravity. Our validity criterion 
can then be summarized as follows: a solution of semi- 
classical gravity is valid when it is stable with respect to 
quantum metric perturbations. This criterion implies to 
consider the quantum correlation functions of the metric 
perturbations. 

It is important to emphasize that the above validity 
criterion incorporates in a unified and self-consistent way 
the two main ingredients of the criteria exposed above. 
Namely, the criteria based on the quantum stress tensor 
fluctuations of the matter fields, and the criteria based 
on the stability of semiclassical solutions against classi- 
cal metric perturbations. In the following discussion we 
will argue that the former is incorporated through the so 
called induced fluctuations and the later though the so 
called intrinsic fluctuations. These correspond to Ford's 
"passive" and "active" fluctuations, respectively. We will 
see that symmetrized quantum two-point metric fluctua- 
tions can always be decomposed as a sum of induced and 
intrinsic fluctuations. 

The paper is organized as follows. In section [H] we 
briefly review the main ingredients of semiclassical grav- 
ity. In section IIIII we introduce stochastic gravity as a 
theory that goes beyond semiclassical theory by incor- 



porating the fluctuations of the quantum stress tensor 
operator. In section HVI our validity criterion is applied 
to the study of flat spacetime as a solution of semiclas- 
sical gravity. The problem of the runaway solutions and 
methods to deal them is discussed. Throughout the pa- 
per in order to emphasize the qualitative aspects we use 
a simplified notation without tensorial indices and for a 
few points we also use qualitative arguments and order of 
magnit ude estimates. We refer the reader to the papers 
|26Ll27ll2a . l29| were the technical details, as well as many 
subtleties that cannot be summarized here, are provided. 
Our metric and curvature conventions are those of Ref. 
|3C|. and we use h = c = 1. 

II. SEMICLASSICAL GRAVITY 

At present semiclassical gravity cannot be rigorously 
derived, but, it can be formally justified in several ways. 
One of them is the leading order in the large N expansion 
|3lj | , where N is the number of independent free quantum 
fields which interact with gravity only. In this limit, after 
path integration one arrives at a theory in which formally 
the gravitational field can be treated as a c-number and 
the quantum fields are fully quantized. 

Semiclassical gravity can be summarized as follows. 
Let g be the metric tensor and (f> a scalar field opera- 
tor. The semiclassical Einstein equation as the dynami- 
cal equation that describes the back-reaction of quantum 
matter on the metric g can be written as 

G a = n{f R ) g , (1) 

where T = T[<p 2 ] is the matter stress tensor in a simpli- 
fied notation, which is quadratic in the field operator </>, 
and k — 8ttG, where G is Newton's constant. This opera- 
tor, being the product of distribution valued operators, is 
ill defined and needs to be regularized and renormalized, 
the R in T R means that the operator has been renormal- 
ized. The angle brackets on the right hand side mean 
that the expectation value of the stress tensor operator 
is computed in some quantum state, say \ip), compatible 
with the geometry described by the metric g. On the 
left hand side G g stands for the Einstein tensor of the 
metric g together with the cosmological constant term 
and other terms quadratic in the curvature which are 
generally needed to renormalize the matter stress tensor 
operator. The quantum field operator <f> propagates in 
the background defined by the metric g, it thus satisfies 
a Klein-Gordon equation, 

(□ 5 -m 2 )0 = O, (2) 

where D g stands for the D'Alambert operator in the 
background of g and m is the mass of the scalar field. 
A solution of semiclassical gravity consists of the set 
(g, (f>, \ip)) where g is a solution of Eq. 0J, cf) is a solu- 
tion of Eq. @ and is the quantum state in which 
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the expectation value of the stress tensor in Eq. (JTJ is 
computed. 

As we recalled in the introduction this theory is in some 
sense unique as a theory that describes the interaction 
of a classical gravitational field with quantum matter. 
As an effective theory it should break down at Planck 
scales. Also, from the right hand side of the semiclassical 
Einstein equation it seems clear that the theory should 
also break down when the fluctuations of the quantum 
stress tensor are large. This has been emphasized by 
Ford and collaborators, and may be illustrated by the 
example of Ref. |(| as follows. 

Let us assume a quantum state formed by an isolated 
system which consists of a superposition with equal am- 
plitude of one configuration with mass Mi and another 
with mass M^. Semiclassical theory as described in Eq. 
(JTJ predicts that the gravitational field of this system 
is produced by the average mass [M\ + Ma)/2, that is 
a test particle will move on the background spacetime 
produced by such a source. However one would expect 
that if we send a succession of test particles to probe the 
gravitational field of the above system half of the time 
they would react to the field of a mass Mi and the other 
half to the field of a mass Mi- If the two masses dif- 
fer substantially the two predictions are clearly different, 
note that the fluctuations in mass of the quantum state 
is of the order of {M\ — M2) 2 . Although the previous 
example is suggestive a word of caution should be said in 
order not to take it too literary. In fact, if the previous 
masses are macroscopic the quantum system decoheres 
very quickly and instead of a pure quantum state it 
is described by a density matrix which diagonalizes in a 
certain pointer basis. Thus for observables associated to 
this pointer basis the matrix density description is equiv- 
alent to that provided by a statistical ensemble. In any 
case, however, from the point of view of the test parti- 
cles the predictions differ from that of the semiclassical 
theory. 

III. STOCHASTIC GRAVITY 

The purpose of stochastic (semiclassical) gravity is to 
be able to deal with the situation of the previous example 
when the predictions of the semiclassical theory may be 
inaccurate. Consequently, our first point is to character- 
ize the quantum fluctuations of the stress tensor. 

The physical observable that measures these fluctua- 
tions is (T ) — (T) 2 . To make this more precise let us 
introduce the tensor operator t = T — (T)I, where I is 
the identity operator, then we introduce the noise ker- 
nel as the four-index bi-tensor defined as the expectation 
value of the anticommutator of the operator t: 

N(x,y) = ±({i(x),i(y)}) g . (3) 

Thus, the noise kernel is the symmetrized connected part 
of the two-point quantum correlation function of the 



stress tensor operator with respect to the state of the 
matter fields. The subindex g here means that this ex- 
pectation value in taken in a background metric g. An 
important property of the symmetric bi-tensor N(x, y) 
is that it is finite because the tensor operator t is finite 
since the ultraviolet divergences of T are cancelled by the 
substraction of (T). Since the operator T is selfadjoint 
N{x, y), which is the expectation value of an anticommu- 
tator, is real and positive semi-definite |28|. Thus, when 
considering the inverse kernel N~ l (x,y), one must work 
in the subspace obtained from the eigenvectors which 
have strictly positive eigenvalues when the noise kernel is 
diagonalized. The last property allows for the introduc- 
tion of a classical Gaussian stochastic tensor £ defined by 

(£(*)>. =0, =N(x,v)- (4) 

This stochastic tensor is symmetric and divergenceless, 
V ■ £ = 0, as a consequence of the fact that the stress 
tensor operator is divergenceless. The subindex s means 
that the expectation value is just a classical stochastic 
average. Note that we assume that £ is Gaussian just 
for simplicity in order to include the main effect of the 
quantum fluctuations. 

The idea now is simple we want to modify the semiclas- 
sical Einstein equation Q by introducing a linear correc- 
tion to the metric tensor g, such as g + h, which accounts 
consistently for the fluctuations of the stress tensor. The 
simplest equation is, 

G g+h = n((f R ) g+h +0, (5) 

where g is assumed to be a solution of equation Q. 
This stochastic equation must be thought of as a lin- 
ear equation for the metric perturbation h which will 
behave, consequently, as a stochastic field tensor. Note 
that the tensor £ is not a dynamical source, since it has 
been defined in the background metric g which is a so- 
lution of the semiclassical equation. Note also that this 
source is divergenceless with respect to the metric, and 
it is thus consistent to write it on the right hand side of 
the Einstein equation. This equation is gauge invariant 
with respect to diffeomorphisms defined by any field on 
the background spacetime [33j. If we take the statistical 
average of equation (J5J it becomes just the semiclassical 
equation for the perturbed metric g + h where now the 
expectation value of T is taken in the perturbed space- 
time. 

The stochastic equation (jSJ) is known as the Einstein- 
Langevin equation. To linear order in h we have [33|, 

(f R ) g+h (x) = -2 J H(x, x') ■ h(x'), (6) 

where the kernel H(x,x') has three terms, one of them 
is proportional to the imaginary part of the expecta- 
tion value of the time ordered two-point stress tensor, 
lm(T(T(x)T(x'))), the second term is proportional to 
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the expectation value of the stress tensor commutator, 
([T(x), T(x')]), and the third is proportional to the func- 
tional derivative of (T) with respect to the metric (ex- 
cluding the implicit dependence on the metric of the field 
4>). Of course, this kernel is also the main ingredient 
of the linearized semiclassical Einstein equation around 
a given background metric g. The other key ingredi- 
ent in the Einstein-Langevin equation is the noise ker- 
nel N(x,y) which defines the stochastic inhomogeneous 
source of the equation. This kernel should be thought 
of as a distribution function, the limit of coincidence 
points has meaning only in the sense of distributions. Ex- 
plicit expressions of this kernel in terms of the two point 
Wightman functions are given in Ref. (3j| on a general 
background. Detailed expressions for this kernel in the 
Minkowski background are given in Ref. [34|], and ex- 
pression based on point-splitting methods have also been 
given in Refs. |9l l35| in other backgrounds. 

The Einstein-Langevin equation has been previously 
derived making use of a formal analogy with open quan- 
tum systems and employing the influence functional for- 
malism |36L l37| . The basis for this approach is a func- 
tional formalism known as closed time path, first intro- 
duced by Schwinger [3^, |3j| |4jj, which is an effective 
action method suitable to derive dynamical equations 
for expectation values of quantum operators; rather than 
transition elements as in the standard effective action 
method. The closed time path formalism was later ap- 
plied to the problem of back-reaction of quantum fields 
on the spacetime metric 0, 0> E3 > m order to derive 
semiclassical Einstein equations. The formalism was then 
applied along the lines of the influence functional for- 
malism to derive Einstein-Langevin equations in several 
contexts Q El |M III El In Ref. [H the 

Einstein-Langevin equation was derived by an axiomatic 
approach by arguing that it is the only consistent gen- 
eralization of the semiclassical Einstein equation which 
takes into account the back-reaction of the matter stress 
tensor fluctuations to lowest order. We have summarized 
the axiomatic approach in this section. 

The solution of the Einstein-Langevin equation JSJ, 
taking into account Eq. |JB}, may be expressed as, 

h(x) = h°(x) + K J G R {x,x')-^x'), (7) 

where hP is a solution of the homogeneous part of equa- 
tion J5J which contains all the information on the initial 
conditions, and G R (x, x') is the retarded propagator with 
vanishing initial conditions associated with the equation 
The two-point correlation function for the metric 
perturbation which is the physically most relevant ob- 
servable can then be written as: 

(h(x)h(y)) s = (h°(x)h°(y)) s + 

k 2 J G R (x, x 1 ) ■ N(x', y') ■ G R (y, y'),(8) 

where the first average, (h°(x)h°(y)) s , is taken with re- 
spect to the initial conditions. 



It turns out that going to leading order in l/N, in the 
large N expansion, one can show that the stochastic cor- 
relation functions for the metric perturbations obtained 
from the Einstein-Langevin equation coincide with the 
symmetrized two-point quantum correlation functions of 
the metric perturbations. The details of the derivation 
will be given in Ref. |5l| and are summarized in Ref. |26| 
for the particular case of a Minkowski background, to 
which we will restrict in section IIVI In this case k in 
Eq. J7J) has to be replaced by the rescaled gravitational 
coupling constant R, — Nk, and the noise kernel for a 
single field N(x,y) must be replaced by (l/N)N(x,y). 
Thus, we have that the symmetrized two-point quantum 
correlation function for the metric perturbation is 

±{{h(x),h(y)}) = {h(x)h(y)) s . (9) 

where the Lorentz gauge condition V ■ {h— (1 / 2)r/Tih) = 
(rj is the Minkowski metric) as well as some initial con- 
dition to fix completely the remaining gauge freedom of 
the initial state should be implicitly understood. 

It should be emphasized that there are two different 
contributions to the symmetrized quantum correlation 
function, which are clearly distinguished in Eq. ||SJ). The 
first contribution is related to the quantum fluctuations 
of the initial state of the metric perturbations and corre- 
sponds to the so called intrinsic fluctuations; here the 
stochastic average must be taken with respect to the 
Wigner distribution function that describes the initial 
quantum state. The second contribution is proportional 
to the noise kernel, it accounts for the fluctuations of the 
stress tensor of the matter fields and corresponds to the 
so called induced fluctuations. These two contributions 
to the two-point correlation functions is also seen in the 
description of some quantum Brownian motion models 
which are typically used as paradigms of open quantum 
systems jU El E3- Both, the intrinsic and induced 
fluctuations, play a role in our stability criterion for the 
solutions of semiclassical gravity. 

The full two-point quantum correlation function for 
the metric (h(x)h(y)} can, in fact, be obtained from the 
Einstein-Langevin equation. Since this correlation can 
be given in terms of the antisymmetrized and the sym- 
metrized quantum correlation function we only need the 
commutator that to leading order in 1 /N is independent 
of the initial state of the metric perturbation and is given 
by 

±({h(x),h(y)})=i K [G R (y,x)-G R (x,y)], (10) 

Note that the information on the retarded propagator is 
already in the linearized semiclassical Einstein equation. 
That is, Eq. © without the stochastic source. 

A. A toy model 

To justify Eq. Q which plays an essential role in our 
criteria for the validity of semiclassical gravity it is useful 
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to introduce a simple toy model for gravity which min- 
imizes the technical complications. The model is also 
useful to clarify the role of the noise kernel and illus- 
trate the relationship between the semiclassical, stochas- 
tic and quantum descriptions. Let us assume that the 
gravitational equations are described by a massless scalar 
field h whose source is another massless scalar field <j) 
which satisfies the Klein-Gordon equation in flat space- 
time □(/> = 0. The field stress tensor is quadratic in the 
field, and independent of h. The classical gravitational 
field equations will be given by 

ah = kT, (11) 

where T is now the (scalar) trace of the stress tensor. 
Note that this is not a self-consistent theory since <f> does 
not react to the gravitational field h. This model obvi- 
ously differs from the standard linearized theory of grav- 
ity discussed previously, where T is also linear in h, but 
it captures some of its key features. 

In the Heisenberg representation the quantum scalar 
field h satisfies 

ah = nf. (12) 

Since T is quadratic in the field operator <f> some regu- 
larization procedure has to be assumed in order for Eq. 
(|12J) to make sense. Since we work in flat spacetime we 
may simply use a normal ordering prescription to reg- 
ularize the operator T. The solutions of this equation, 
i.e. the field operator at the point x, which we call h x 
in this subsection to avoid confusion with the more stan- 
dard notation, h(x), used in the rest of the paper, may 
be written in terms of the retarded propagator G xx i of 
the D'Alambertian as, 

h x = h°+K J G xx -Tx>, (13) 

where h x is the free field which carries information on 
the initial conditions and the state of the field. From 
this solution we may compute, for instance, the symmet- 
ric two-point quantum correlation function (the anticom- 
mutator) 

({h x ,hy}) = ({h° x ,h° y })+K 2 j G xx ,G yy ,({f x ,,fy,}), 

where the expectation value is taken with respect to the 
quantum state in which both fields <f> and h are quantized. 
We have assumed (h°) = for the free field. 

We can now consider the semiclassical theory for this 
problem. If we assume that h is classical and the mat- 
ter field is quantum the semiclassical limit may just be 
described by substituting into the classical equation i|ll|) 
the stress trace by the expectation value of the scalar 
stress operator (T) , in some quantum state of the field <fi. 
We may simply renormalize the expectation value of T 
using normal ordering, then for the vacuum state of the 



field <f>, we would simply have (T% = 0. The semiclassical 
theory thus reduces to 

□/i = k(T). (15) 

The two point function h x h y that one may derive from 
this equation depends on the two point function (T x )(T y ) 
and clearly cannot reproduce the quantum result of 
Eq. I|14f> which depends on the expectation value of two- 
point operator ({T x ,T y }). That is, the semiclassical the- 
ory entirely misses the fluctuations of the scalar stress 
operator T. 

To extend this semiclassical theory in order to account 
for such fluctuations, we introduce the noise kernel as we 
did in the previous section. Thus, we define 

N xy = ~({i x j y }) (16) 

where t = T — (T), and we have used again the sub- 
index notation to avoid confusion with the noise kernel 
of the previous section. The bi-scalar N xy is real and 
positive-semidcfinite, as a consequence of t being self- 
adjoint |28j . Consequently we can introduce a Gaussian 
stochastic field as: 

(0s=0, (Uv). = N xv . (17) 

where the subscript s means a statistical average. 

The extension of the semiclassical equation may be 
simply performed by adding to the right-hand side of 
the semiclassical equation (|15|) the stochastic source £, 
which accounts for the fluctuations of T as follows, 

Oh = K ((T ) + £) . (18) 

This equation is in the form of a Langevin equation: the 
field h is classical but stochastic and the observables we 
may obtain from it are correlation functions for h. In fact, 
the solution of this equation may be written in terms of 
the retarded propagator as, 

h x = h° x + K J G xx , (<T X /) + iJ) , (19) 

from where the two point correlation function for the 
classical field h, after using the definition of £ and that 
(h°) s = 0, is given by 

{hxhy). = (h° x h° y ) s + y J G xx ,G yy ,({f x ,,f y ,}). (20) 

Note that in writing (••■)„ here we are assuming a double 
stochastic average, one is related to the stochastic process 
£ and the other is related to the free field h° which is 
assumed also to be stochastic with an initial distribution 
function to be specified. 

Comparing Eqs. (|14|l and (121 III we see that the respec- 
tive second term on the right-hand side are identical (ex- 
cept for a factor of 2 due to the symmetrization) provided 
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the expectation values are computed in the same quan- 
tum state for the field <j>. The fact that the field h is also 
quantized in (|14|l does not change the previous state- 
ment; recall that T does not depend on h. The nature of 
the first term on the right-hand sides of equations 114|) 
and l|20|) is different: in the first case it is the two-point 
quantum expectation value of the free quantum field h 
whereas in the second case it is the stochastic average of 
the two point classical homogeneous field h°, which de- 
pends on the initial conditions. Now we can still make 
these terms equal to each other (with the factor of 2) 
if we assume for the homogeneous field h° a Gaussian 
distribution of initial conditions such that 

(KK)s = \{{hl,hl)). (21) 

This Gaussian stochastic field h° can always be defined 
due to the semi-positivity of the anti-commutator. Thus, 
under this assumption on the initial conditions for the 
field h the two point correlation function of Eq. I|20(l 
equals the quantum expectation value of Eq. (|14|l exactly. 
Thus, we have 

\{{h x ,h y }) = {h x h y ) s , (22) 

which may be compared to Eq. 10. Comparing with 
the linearized theory of gravity described in the previous 
section we see that (T) depends also on h, both explic- 
itly and also implicitly through the coupling of <f> with 
h. The retarded propagator here G xx > is then replaced 
by the propagator Gn(x,x r ) of the previous section and 
the functions which are here the free metric pertur- 
bations are replaced by the homogeneous solutions of the 
previous section. 



IV. STABILITY OF FLAT SPACETIME 

Let us now apply our validity criterion to flat space- 
time. One particularly simple and interesting solution of 
semiclassical gravity is the Minkowski metric. In fact, 
when the quantum fields are in the Minkowski vacuum 
state one may take the renormalized expectation value 
of the stress tensor (T R ) = (this is equivalent to as- 
suming that the cosmological constant is zero) and the 
Minkowski metric 77 is a solution of the semiclassical Ein- 
stein equation Thus, we can look for the stability of 
flat spacetime against quantum matter fields. According 
to the criteria we have established we have to look for the 
behavior of the two-point quantum correlations for the 
metric perturbations h over the Minkowski background 
which are given by Eqs. © and 10. As we have em- 
phasized several times these fluctuations separate in two 
parts: the first term on the right hand side of Eq. (JHJ) 
corresponds to the intrinsic fluctuations, and the second 
term corresponds to the induced fluctuations. 



A. Intrinsic fluctuations 

Let us first consider the intrinsic fluctuations, 

(h°(x)h°(y)) s , (23) 

where h° are the homogeneous solutions of the Einstein- 
Langevin equation (jSJ), or equivalently the linearly per- 
turbed semiclassical equation, and the statistical average 
is taken with respect to the Wigner distribution that de- 
scribes the initial quantum state of the metric pertur- 
bations. Since these solutions are described by the lin- 
earized semiclassical equation around flat spacetime we 
can make use of the results derived in Refs. [Tit EH l23j| . 
The solutions for the case of a massless scalar field were 
first discussed in Ref. and an exhaustive description 
can be found in Appendix A of Ref. [2l| • Decomposing 
the metric perturbation into scalar, vectorial and tenso- 
rial parts and computing the linearized Einstein tensor, 
one gets a vanishing result for the vectorial part of the 
metric perturbation; the scalar and tensorial components 
of the metric perturbation give rise, respectively, to the 
scalar and tensorial components of the linearized Einstein 
tensor. The vectorial part is found to vanish whereas the 
scalar and tensorial contributions for a massless and con- 
formally coupled scalar field (see Ref. [2lJ for the mass- 
less case with arbitrary coupling and Refs. |U Q for 
the general massive case) satisfy the following equations: 

(l + 12K/3p 2 )G (s) (p) =0, (24) 

J&O + iap*?)*^- - (25) 

where in the last equation the prescription that the time 
component of p has a small imaginary part, p° + ie, is 
taken. Here G{p) stands for the Fourier transform of 
the linearized Einstein tensor, the upper indices S and T 
stand for scalar and tensorial respectively, [3 is a dimcn- 
sionless renormalized parameter that multiplies some of 
the quadratic terms in the curvature in the effective ac- 
tion for the gravitational field, and /z is a renormalization 
mass scale. See Ref. 26] for a more complete description. 

For the scalar component when (3 = the only solution 
is & s \p) = 0. When j3 > the solutions for the scalar 
component exhibit an oscillatory behavior in spacetime 
coordinates which corresponds to a massive scalar field 
with to 2 = (12k|/3|) _1 ; for (3 < the solutions corre- 
spond to a tachyonic field with m? = — (12/c|/3|) : in 
spacetime coordinates they exhibit an exponential behav- 
ior in time, growing or decreasing, for wavelengths larger 
than 47t(3k|/?|) 1 / 2 , and an oscillatory behavior for wave- 
lengths smaller than 47t(3k|/3|) 1 / 2 . On the other hand, 
the solution G^^p) = is completely trivial since any 
scalar metric perturbation h(p) giving rise to a vanishing 
linearized Einstein tensor can be eliminated by a gauge 
transformation as explained in Ref. p6| . 
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For the tensorial component, when /i < /i cr it = 
«- 1 (1207r) 1 / 2 e'>' ) where l p is the Planck length (I* = k/8tt) 
the first factor in Eq. Q25J1 vanishes for four complex val- 
ues of p° of the form ±uj and ±uj*, where is some 
complex value. We will consider here the case in which 
M < Merit; a detailed description of the situation for 
M > Merit can be found in Appendix A of Ref . [2lJ . The 
two zeros on the upper half of the complex plane cor- 
respond to solutions in spacetime coordinates exponen- 
tially growing in time, whereas the two on the lower half 
correspond to solutions exponentially decreasing in time. 
Strictly speaking, these solutions only exist in spacetime 
coordinates, since their Fourier transform is not well de- 
fined. They are commonly referred to as runaway solu- 
tions and for \i ~ l^ 1 they grow exponentially in time 
scales comparable to the Planck time. 

In order to deal with those unstable solutions, one pos- 
sibility is to employ the order reduction prescription [22j , 
which we will briefly summarize in the last subsection. 
With such a prescription we are left only with the so- 
lutions which satisfy G(p) = 0. The solutions for h(p) 
simply correspond to free linear gravitational waves prop- 
agating in Minkowski spacetime expressed in the trans- 
verse and traceless (TT) gauge. When substituting back 
into Eq. ill'-'il) and averaging over the initial conditions we 
simply get the symmetrized quantum correlation func- 
tion for free gravitons in the TT gauge for the state given 
by the Wigner distribution. As far as the intrinsic fluc- 
tuations are concerned, it seems that the order reduc- 
tion prescription is too drastic, at least in the case of 
Minkowski spacetime, since no effects due to the interac- 
tion with the quantum matter fields are left. 

A second possibility, proposed by Hawking et al. 
[HBl |5(| , is to impose boundary conditions which discard 
the runaway solutions that grow unbounded in time and 
correspond to a special prescription for the integration 
contour when Fourier transforming back to spacetime co- 
ordinates. Following that procedure we get, for example, 
that for a massless conformally coupled scalar field with 
[3 > the intrinsic contribution to the symmetrized quan- 
tum correlation function coincides with that of free gravi- 
tons plus an extra contribution for the scalar part of the 
metric perturbations which renders Minkowski spacetime 
stable but plays a crucial role in providing a graceful exit 
for inflationary models driven by the vacuum polariza- 
tion of a large number of conformal fields. Such a massive 
scalar field would not be in conflict with present obser- 
vations because, for the range of parameters considered, 
the mass would be far too large to have observational 
consequences |55j . 



B. Induced fluctuations 

The induced fluctuations are described by the second 
term in Eq. ©. They are induced for the noise kernel 
that describes the stress tensor fluctuations of the matter 



fields, 

^ J G R (x, x') ■ N(x', y') ■ G R {y, y% (26) 

where we write the expression in the large N limit. The 
contribution corresponding to the induced quantum fluc- 
tuations is equivalent to the stochastic correlation func- 
tion obtained by considering just the inhomogeneous part 
of the solution to the Einstein-Langevin equation: the 
second term on the right-hand side of Eq. © . Taking all 
that into account, it is clear that we can make use of the 
results for the metric correlations obtained in Ref. [34[ 
by solving the Einstein-Langevin equation. In fact, one 
should simply take N = 1 to transform our expressions 
to those of Ref. [34| and, similarly, multiply the noise 
kernel in the expressions of that reference by N so that 
they can be used here, which follows from the fact that 
we have ./V independent matter fields. 

The same kind of exponential instabilities in the run- 
away solutions of the homogeneous part of the Einstein- 
Langevin equation also arise when computing the re- 
tarded propagator Gr . In order to deal with those insta- 
bilities, similar to the case of the intrinsic fluctuations, 
one possibility is to make use of the order reduction pre- 
scription. The Einstein-Langevin equation becomes then 
G(p) = R,£,(p). The second possibility, following the pro- 
posal of Hawking et al., is to impose boundary condi- 
tions which discard the exponentially growing solutions 
and translate into a special choice of the integration con- 
tour when Fourier transforming back to spacetime co- 
ordinates the expression for the propagator. In fact, it 
turns out that the propagator which results from adopt- 
ing that prescription coincides with the propagator that 
was employed in Ref. |34j . Note, however, that this prop- 
agator is no longer a strictly retarded propagator since 
it exhibits causality violations at Planck scales. A more 
detailed discussion on all these points can be found in 
Appendix E of Re f. El . 

Following Ref. |34j . the Einstein-Langevin equation 
can be entirely written in terms of the linearized Ein- 
stein tensor. The equation involves second derivatives of 
that tensor, and in terms of its Fourier components G(p) 
takes the form 

(l + F(p)).G(p) = «|(p), (27) 

where F is a four- index tensor which depends on p 2 \np 2 
when the field is massless and conformally coupled. This 
reflects the fact that we have second derivatives of the 
Einstein tensor and the nonlocality of the Einstein- 
Langevin equation (or also of the perturbed semiclassical 
equation) . From equation l|27|l one may obtain the corre- 
lation functions for G(p), (G(p)G(q)) s , which are invari- 
ant under gauge transformations of the metric pertur- 
bations. Writing the linearized Einstein tensor in terms 
of the metric perturbation, which takes a particularly 
simple form in the Lorentz gauge, one may derive the 
correlation functions for h(p): (h(p)h(q)) s . Finally, the 
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correlation functions in spacetime coordinates can be eas- 
ily obtained by Fourier transforming these correlations. 
For massless and conformally coupled matter fields ex- 
plicit results are given in Ref. [34| , they have the general 
expression: 

R 2 f e lp -( x -y)pB(— p 2 ) 

(h(x)h(y)) s = 72Q ^ N j | 1 + ( S / 9607r 2) p2 l n ( p 2 ^2)12 

(28) 

where P is a four-index projection tensor. This corre- 
lation function for the metric perturbations is in agree- 
ment with the real part of the propagator obtained by 
Tomboulis in Ref. |53 using a large N expansion. 

To estimate this integral let us consider spacelike sep- 
arated points x — y — (0, r) and introduce the Planck 
length I p . It is not difficult to see that for space 
separations |r| » l p we have 

(h(x)h(y)) s ~ £ (29) 

and for |r| ~ Nl p we have 

(h(x)h(y)) s ~ e-W/'p-^. (30) 

Since these fluctuations are induced by the matter stress 
fluctuations we infer that the effect of the matter fields 
is to suppress metric fluctuations at small scales. On the 
other hand, at large scales the induced metric fluctua- 
tions are small compared to the free graviton propagator 
which goes like ip/|r| 2 . 

We thus conclude that, once the instabilities giving rise 
to the unphysical runaway solutions have been discarded, 
the fluctuations of the metric perturbations around the 
Minkowski spacetime induced by the interaction with 
quantum scalar fields are indeed stable. Instabilities 
would lead to a divergent result when Fourier transform- 
ing back to spacetime coordinates. Note that when the 
order reduction prescription is used the p 2 hip 2 terms are 
absent in the corresponding Eq. (|28H . Thus, in contrast 
to the intrinsic fluctuations, there is still a nontrivial con- 
tribution to the induced fluctuations due to the quantum 
matter fields in this case. 



C. Order reduction prescription and large N 

Runaway solutions are a typical feature of equations 
describing back-reaction effects, such is in classical elec- 
trodynamics, and are due to higher than two time deriva- 
tives in the dynamical equations. In a very schematic way 
the semiclassical Einstein equations have the form 

G h + %G h = 0, (31) 

where Gh stands for the linearized Einstein tensor over 
the Minkowski background, say, and we have simplified 
the equation as much as possible. The second term of 



the equation is due to the vacuum polarization of mat- 
ter fields and contains four time derivatives of the metric 
perturbation. Some specific examples of such an equation 
are, in momentum space, Eqs. (|24|l and l|25|l . The order 
reduction procedure is based on treating perturbatively 
the terms involving higher order derivatives, differenti- 
ating the equation under consideration and substituting 
back the higher derivative terms in the original equa- 
tion keeping only terms up to the required order in the 
perturbative parameter. In the case of the semiclassi- 
cal Einstein equation, the perturbative parameter is I 2 . 
If we differentiate twice Eq. 131|) with respect to time 
it is clear that the second order derivatives of the Ein- 
stein tensor are of order I 2 . Substituting back into the 
original equation, we get the following equation up to 
order l p : Gh = + 0(1*). Now, there are certainly no 
runaway solutions but also no effect due to the vacuum 
polarization of matter fields. Note that the result is not 
so trivial when there is an inhomogeneous term on the 
right hand side of Eq. I|31|) , this is what happens with the 
induced fluctuations predicted by the Einstein-Langevin 
equation. 

Semiclassical gravity is expected to provide reliable re- 
sults as long as the characteristic length scales under con- 
sideration, say L, satisfy that L ^> l p |2l|. This can 
be qualitatively argued by estimating the magnitude of 
the different contributions to the effective action for the 
gravitational field, considering the relevant Feynman dia- 
grams and using dimensional arguments. Let us write the 
effective gravitational action, again in a very schematic 
way, as 

Ses = J V=g (^R + olR 2 + i 2 p R 3 + ...), (32) 

where R is the Ricci scalar. The first term is the usual 
classical Einstein-Hilbert term, the second stands for 
terms quadratic in the curvature (square of Ricci and 
Weyl tensors) this terms appear as radiative corrections 
due to vacuum polarization of matter fields, here a is 
an dimensionless parameter presumably of order 1, the 
R 3 terms are higher order corrections which appear for 
instance when one considers internal graviton propaga- 
tors inside matter loops. Let us assume that R ~ L~ 2 
then the different terms in the action are of the order 
of R 2 ~ L~ 4 and l 2 R 3 ~ lpL~ 6 . Consequently when 
L ^> 1^, the term due to matter loops is a small correc- 
tion to the Einstein-Hilbert term (l/l 2 )R ^> R 2 , and this 
term can be treated as a perturbation. The justification 
of the order reduction prescription is actually based on 
this fact. Therefore, significant effects from the vacuum 
polarization of the matter fields are only expected when 
their small corrections accumulate in time, as would be 
the case, for instance, for an evaporating macroscopic 
black hole all the way before reaching Planckian scales. 

However if we have a large number N of matter fields 
the estimates for the different terms change in a remark- 
able way. This is interesting because the large N expan- 
sion seems the best justification for semiclassical grav- 
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ity. In fact, now the vacuum polarization terms involv- 
ing loops of matter are of order NR 2 ~ NL~ A . For this 
reason the contribution of the graviton loops, which is of 
order R 2 , can be neglected in front of the matter loops; 
this justifies the semiclassical limit. Similarly higher or- 
der corrections are of order Nl 2 R 3 ~ Nl 2 L~ e . Now there 

is a regime, when L ~ y/~Nl p , where the Einstein-Hilbcrt 
term is comparable to the vacuum polarization of mat- 
ter fields, (l/l 2 )R ~ NR 2 , and yet the higher correction 
terms can be neglected because we still have L ^> l p , pro- 
vided N ^> 1. This is the kind of situation considered 
in trace anomaly driven inflationary models 1551 . such 
as that originally proposed by Starobinsky |58j |. where 
the exponential inflation is driven by a large number of 
massless conformal fields. The order reduction prescrip- 
tion would completely discard the effect from the vacuum 
polarization of the matter fields even though it is compa- 
rable to the Einstein-Hilbert term. In contrast, the pro- 
cedure proposed by Hawking et al. keeps the contribution 
from the matter fields. Note that here the actual physi- 
cal Planck length l p is considered, not the rescaled one, 
lp = k/8tt, which is related to l p by I 2 — k/8it — l 2 /N. 

V. CONCLUSIONS 

An analysis of the stability of any solution of semiclas- 
sical gravity with respect to small quantum corrections 
should consider not only the evolution of the expectation 
value of the metric perturbations around that solution, 
but also their fluctuations, encoded in the quantum corre- 
lation functions. Making use of the equivalence (to lead- 
ing order in 1/N, where N is the number of matter fields) 
between the stochastic correlation functions obtained in 
stochastic semiclassical gravity and the quantum correla- 
tion functions for metric perturbations around a solution 
of semiclassical gravity, the symmetrized two-point quan- 
tum correlation function for the metric perturbations can 
be decomposed into two different parts: the intrinsic fluc- 
tuations due to the fluctuations of the initial state of 



the metric perturbations itself, and the fluctuations in- 
duced by their interaction with the matter fields. If one 
considers the linearized perturbations of the semiclassi- 
cal Einstein equation, only information on the intrinsic 
fluctuations can be retrieved. On the other hand, the in- 
formation on the induced fluctuations naturally follows 
from the solutions of the Einstcin-Langevin equation. 

As a specific example, we have analyzed the sym- 
metrized two-point quantum correlation function for the 
metric perturbations around the Minkowski spacetime in- 
teracting with N scalar fields initially in the Minkowski 
vacuum state. Once the ultraviolet instabilities which are 
ubiquitous in semiclassical gravity and are commonly re- 
garded as unphysical, have been properly dealt with by 
using the order reduction prescription or the procedure 
proposed by Hawking et al. ,55, ' 56} , both the intrinsic 
and the induced contributions to the quantum correla- 
tion function for the metric perturbations are found to 
be stable (2||. Thus, we conclude that Minkowski space- 
time is a valid solution of semiclassical gravity. 
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